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SUMMARY 


This Final Report represents the final contract deliverable for NASA MSFC Purchase Order No. 
H-13047D, entitled: 

Design Study - Rocket Based MHD Generator 

The Final Report addresses the technical feasibility and design of a rocket based MHD generator 
using a sub-scale LOx/RP rocket motor. The design study was constrained by assuming the 
generator must function within the performance and structural limits of an existing magnet and 
by assuming realistic limits on (1) the axial electric field, (2) the Hall parameter, and (3) current 
density, and (4) heat flux (given the criteria of heat sink operation). The major results of the 
work are summarized as follows: 

• A Faraday type of generator with rectangular cross section is designed to operate with a 
combustor pressure of 300 psi. Based on a magnetic field strength of 1.5 Tesla, the 
electrical power output from this generator is estimated to be 54.2 KW with potassium 
seed (weight fraction 3.74%) and 92 KW with cesium seed (weight fraction 9.66%). 
The former corresponds to a enthalpy extraction ratio of 2.36% while that for the latter 
is 4. 16%. 

• A conceptual design of the Faraday MHD channel is proposed, based on a maximum 
operating time of 10 to 15 seconds. This concept utilizes a phenolic back wall for 
inserting the electrodes and inter-electrode insulators. Copper electrode and aluminum 
oxide insulator are suggested for this channel. 

• A testing configuration for the sub-scale rocket based MHD system is proposed. 

An estimate of performance of an ideal rocket based MHD accelerator is performed. With a 
current density constraint of 5 Amps/cm 2 and a conductivity of 30 Siemens/m, the push power 
density can be 250, 431, and 750 MW/m 3 when the induced voltage uB have values of 5, 10, and 
15 KV/m, respectively. 
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1. Introduction 


The concept of coupling a rocket motor with magnetohydrodynamic (MHD) devices was 
proposed several decades ago. 1 ' 2 ' 3 ' 4 The is due to not only that rocket motors are an appropriate 
plasma source (after seeding) for MHD applications but also a simple structured, light weight 
MHD device is feasible. Recently, there has been a resurgent interest in MHD aerospace 
applications, in which system weight is always a major concern. Rocket based MHD system is 
considered a good candidate to fit in these aerospace applications. A couple of rocket based 
MHD systems have been proposed. Among them, some were for electric power generation, 
others for flow acceleration, or even a combination of both. 10 

The present rocket based MHD generator design study is deemed as the first phase of a research 
project for developing such a system. Next phase will be fabrication and laboratory test for the 
rocket MHD system design concluded in the present study. Sub-scale rocket motors with 300 or 
1000 psi combustion pressure of NASA Marshall Space Flight Center (MSFC) and a magnet 
system used in early MHD generator studies at Energy Conversion Facility (ECF) of The 
University of Tennessee Space Institute (UTSI) 11 are proposed to be used in the next phase 
laboratory set up. Due to certain practical limitations associated with these rocket motors and 
the magnet system, the rocket based MHD generator designed herein might not show great 
performance. The major goals of the entire research project are to demonstrate power 
production from the designed generator, to set up a data base for verification of analytical 
modeling work, and to accumulate experience in designing and operating an MHD generator of 
this kind. Once these goals are reached, large scale rocket based MHD generators can be later 
designed to provide great performance with confidence. 

This report overviews work performed on the rocket based MHD generator design study 
contract awarded to ERC. Inc. from NASA Marshall Space Flight Center. Possibilities to 
examine rocket based MHD accelerator are resumed through the entire research project. 


2. LOx/RP Rocket Motor Thermo/Plasmadynamic Characterization 

The MHD generator envisioned in this study is based on sub-scale LOx/RP rocket motors 
available at the Marshall Space Flight Center. With the same nozzle geometry, these rocket 
motors can be divided into two categories according to combustion pressures, e.g. 300 and 1000 
psia. Mass flow rates, thrust, exit temperature, and pressure of these two kinds of rocket motors 
burning with 50/50 of LOx/RP were calculated using NASA thermo-chemical equilibrium code 
(SP-273) 12 and listed in Table 1 along with important geometrical factors. The thrust thus 
obtained assumed ideal expansion to vacuum. Sketch of the 500 pound thrust liquid fuel rocket 
motor is shown in Figure 1. 
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Figure 1 Sub-Seale Rocket Motor of MSFC (SNTL - D - 500 - 1000, Liquid Fuel) 



Stagnation Pressure (psi) 

300 

1000 

Stagnation Temperature (K) 

3473 

3634 

Throat Diameter (in) 

0.6 

0.6 

Nozzle Exit Diameter (in) 

1.88 

1.88 

Exit Temperature(K) 

2160 

2126 

Exit Pressure (atm) 

0.295 

0.953 

Exit Mach No. 

3.09 

3.14 

Throat Sonic speed (m/s) 

1170. 

1191 

Mass Row Rate(kg/s) 

0.211 

0.695 

Thrust (lb f ) 

150 

500 


Table 1. Geometry and Theoretical Performance of MSFC’s Sub-scale Rocket Motors 


The primary objective of the present design study is to reach a conceptual design of a rocket 
based MHD system which produces maximum power output based on MSFC’s sub-scale rocket 
motors and ECF magnet system. Treating the MHD generator as additional payload to the 
rocket system, a compact MHD device is preferred. To minimize the weight of an MHD 
channel, the power density, defined as the electrical power output per unit of internal volume, of 
the generator has to be maximized. 

It is widely known that power density of an MHD generator is proportional to the product of the 
plasma conductivity with the squares of the plasma velocity and imposed magnetic field 
strength. That is 

Power Density ~ <7 u B 2 (1) 

where cr is the electrical conductivity of plasma, u is the velocity, and B is the magnetic field 
strength. To maximize MHD performance, the highest power density possible needs to be 
maintained in the channel of an MHD generator. 


2.1 Seeding 

From the temperature range shown in Table 1, a LOx/RP rocket motor requires seeding its 
combustion gas with an easily ionized alkali salt to provide high conductivity plasma. 
Applications of cesium hydroxide (CsOH) or potassium hydroxide (KOH) as the seeding 
component is examined in this design study. Figure 2 presents calculated plasma conductivities 
at the combustor conditions for various weight percentages of potassium and cesium. As 
indicated in the figure, peak value of conductivity exists in each curve. For example, the 
conductivity maximizes at 5.2% and 16.5% for potassium and cesium respectively at combustor 
pressure of 300 psi. Cesium seeding provides higher plasma conductivity than potassium by 
virtue of its lower ionization potential. 
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Figure 2 Plasma Conductivity Variation with Seeding Percentage of Potassium and Cesium at 
Various Combustion Conditions of LOx/RP Rocket Motors. 


On the other hand, cesium is extremely toxic and causes concerns of environmental pollution 
when exhausted by a rocket. Also cesium is considerably more expensive than potassium. 
Therefore, potassium is selected as the primary seeding component herein. Cases using cesium 
as seed would also be considered in this investigation as a maximum effort to achieve a high 
power density by maximizing the plasma conductivity when only limited magnetic field strength 
is available. 


2.2 Generator Entry Point Selection 

It is convenient to further subdivide the expression (1) of power density into the terms B 2 and the 

2 

product — where the latter is termed as the power density parameter. The power density 

4 

parameter is a plasma characteristic. From its definition, the plasma conductivity is an 
exponentially increasing function with temperature while plasma velocity is approximately a 
decreasing square root function with temperature in a nozzle flow. 


4 


ERC-R-97-017 



The potential to produce electrical power in the expanded plasma states is of primary importance 
for a rocket based MHD power system. During the expansion, the plasma temperature and 
pressure drop while its velocity increases. While the plasma conductivity increases with 
decreasing pressure at constant temperature, the temperature dependence is strong and dominant 
so that the conductivity in general decreases as the plasma expands to higher speed. That is, the 
power density parameter is closely related with plasma temperature and exhibits an optimal 
value when plasma flowing through a rocket nozzle. Along the expansion section of nozzle, the 
point where the optimal power density exists is a good choice as the entry of MHD channel to 
ensure maximum MHD performance. 

Figures 3 and 4 are the power density parameter mappings for potassium and cesium 
respectively produced by varying both the seeding percentage and the expansion ratio when the 
combustion pressure is 300 psi.. For potassium, the peak value of power density parameter is 
found at 3.74% seed and expansion ratio of 8.0 while, for cesium, it is located at 9.66% seed and 
expansion ratio of 14.1. The former expansion ratio (for potassium) is less than that of MSFC’s 
rocket but the latter (for cesium) is greater, according to the geometry listed in Table 1. In 
practice, a relatively wide range of expansion ratio can be selected without too much reduction. 
For example, along the curve of 3.74% potassium seeding in Figure 3, the expansion ratio can be 
chosen anywhere between 6 to 10 with less than 1.2% decrease in power density parameter from 
its peak value. From the 9.66% cesium seeding line in Figure 4, any expansion ratio between 
9.5 to 20 can be used if 1.2% reduction in power density parameter relative to its optimal value 
is allowable. 

For 1000 psi combustor pressure, the peak power density parameter occurs at about the same 
seed percentage and expansion ratio as that for 300 psi case for both potassium and cesium. 
Keeping the same seeding percentage, variations in the expansion ratio between 6 to 10 for 
potassium and between 9.5 to 20 for cesium result in less than 1.5% drop in power density 
parameter with respect to its peak value. 


3. Design Analysis 

Two major types of MHD generator in terms of loading configurations are Faraday and Hall 
types, as shown in Figure 5. Electrode pair of a Hall generator is in the flow direction (i.e., x- 
direction) while that of a Faraday generator is perpendicular to the flow train (i.e., y-direction). 
With segmented electrodes, axial electrical currents (also named as Hall currents) inside a 
Faraday generator can ideally be neutralized and only Faraday currents flow through the 
electrodes. With an externally applied magnetic field in the z-direction, the induced Lorentz 
force from MHD interaction will lie in direction against flow train in a Faraday type of generator 
to decelerate the flow and convert the flow kinetic energy into the electrical energy. 
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Parameter (Mw/m3* T 2) | | Power Density Parameter (MW/m3“T2) 



Power Density Parameter Distribution along 300 psi LOx/RP Rocket Nozzle with 
Different Potassium Seeding Percentages. 



Figure 4 Power Density Parameter Distribution along 300 psi LOx/RP Rocket Nozzle with 
Different Cesium Seeding Percentages. 
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Figure 5 Loading Configuration of a Faraday Generator with Segmented Electrodes and a Hall 
Generator 


On the other hand, the wall of a Hall generator is shorted, currents in the axial direction are that 
flow through the electrodes. With a magnetic field in the z-direction, the Lorentz force in a Hall 
generator has a component perpendicular to the flow train and creates a side force against the 
channel wall. This side force, if significant in magnitude, might perturb the flow train, cause 
boundary layer separation, and make a rocket to turn, when a rocket based MHD generator 
employs a Hall configuration. Based on the operating conditions of MSFC’s sub-scale rocket 
motor, an order of magnitude estimate of Lorentz force in a Hall generator is carried out below. 

Axial current density, j x , in a Hall generator is calculated as 

j K =-^—r{\-K H )auB (2) 

1 + Cl>T" 

where an is Hall parameter defined as the ratio between electron cyclotron frequency and 
electron coilisional frequency with atoms, molecules, and ions. Kh is the loading factor of a Hall 
generator. 


To balance the Lorentz force, a pressure force, Ap y , in the y-direction will form as 


Av 


- j, B 


Combining expressions (2) and (3), the order of Ap y can be estimated by 

ClYT 

An « — ( 1 — K„) GuB 1 Av 

1 + an~ 


(3) 

(4) 
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In a sub-scale rocket based MHD generator, typical values for variables in the right hand side of 
expression (4) are cor ~ 5, K H ~ 0.5, <7-20 Siemens/m, u ~ 2500 m/s, B 2 T, and Ay ~ 0.1 m. 
From these, we can get Ap y = 2000 pa. With a density p = 0.03 kg/m 3 , ratio between this 
Lorentz induced pressure force and dynamic pressure in the flow train is about lxlO' 2 . 
Therefore, we can conclude that the side force induced by MHD interaction in a sub-scale rocket 
based Hall generator does not perturb the flow train significantly. 

An MHD generator of Hall configuration has simpler structure compared to one of the Faraday 
type since only one single electrode pair and single load is needed in a Hall generator. This 
simplicity makes Hall configuration more attractive in applications where system weight is of 
primary importance, even though a linear Hall generator is usually less efficient than one of 
Faraday type. On the other hand, it is noticed that the side force induced by Hall current is 
negligible only for small scale MHD generators. For a full scale MHD scale, Ay can be one 
order of magnitude larger and, under such circumstances, the Lorentz induced side force is no 
longer negligible. Therefore, a Faraday type of generator is suggested in the present study. 


3.1 Magnet Constraints 

For cost savings, an existing magnet will be used in the experimental hardware set up following 
this design study. Therefore, it is necessary to size the experiment before detailed calculations 
for the MHD generator channel can be performed. This is contrast to conventional MHD design 
practice. In a conventional design process, the guiding criteria for the generator channel are 
couched in terms of plasma mass flow and electrical conductivity available as well as power 
output sought. Subsequent for the final channel configuration, the magnet is constructed to 
match the required field strength and interaction length. In the present study, the procedure is to 
accommodate an MHD channel within an existing magnetic field envelop. 

The magnet proposed to be used in the follow up experiment is the existing ECF’s magnet at 
UTSI. This magnet is a conventional water-cooled electromagnet with a C-frame core. The 
pole gap is 10 cm and the pole faces were originally rectangular shape (91.44 cm x 15.2 cm) and 
later on modified to a trapezoidal shape. This magnet was originally powered by DC welders 
and provided a maximum field strength of 2.2 Tesla. 

Figure 6 shows the normalized magnetic induction along the magnet axis. The magnetic field 
distribution is relatively flat nearby its peak for a distance approximately 20 to 24 inches. In 
order to maximize power density, which is proportional to the square of the field strength, it is 
desirable to have a uniform high strength field. The generator must be able to fit within the 
volume with uniform field. In the present design, the constant field strength is selected to be 1.5 
Tesla and the effective length of the MHD channel is 20 inches, evenly divided into 20 electrode 
pitches. This level of magnetic field strength is believed to be achievable after rebuilding the 
magnet system, even though it has been out of service since early 80’s. 
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Figure 6 Normalized Magnetic Flux Density Distribution of ECF Magnet 


For sake of machining simplicity and cost saving, an MHD channel with circular cross section 
was thought to be a possible candidate for designing a rocket based MHD generator where 
rocket combustor and nozzle are of circular shape. Reference ? is an example for applying such 
a circular cross section channel. In the present design, due to the narrowness of the magnet pole 
gap (— 4 inches), circular cross section is impossible to feed the necessary mass flow rate, taking 
into considerations of the thickness of channel walls. The generator designed can only have a 
rectangular channel with a constant width but diverging in its height. Matching a cylindrical 
rocket to a rectangular channel requires a transition nozzle which accomplishes the cross- 
sectional shape transition. This transition nozzle has to be carefully designed and fabricated to 
avoid comer interference and subsequently flow non-uniformity. 


3.2 Other Constraints 

For a Faraday generator with segmented electrodes, the electrical field in the axial direction has 
to be limited to avoid happening of electrode arcing phenomena. Once arcing occurs, the 
adjacent electrodes can be shorted to degrade the generator performance. The locally high 
dissipation from highly concentrated arcing discharges might result in erosion of the electrode 
surface and shorten the channel operating life. In the design, the axial electrical field is limited 
to below 40 volts/cm by adjusting the generator loading. 
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The plasma in the MHD channel is better to operate under a critical Hall parameter. If Hall 
parameter exceeds to its critical value, eddy currents appear in plasma and large scale plasma 
non-uniformities form in the channel. These non-uniformities increase the plasma internal 
impedance and subsequently lower the generator performance. The critical Hall parameter is 
selected to be 4. 

Faraday currents generated in the generator flow to the external loads through electrodes. The 
Faraday current density has to be constrained below a certain value, depending on electrode 
material and geometry, to avoid excessive heat generated on electrode surface and subsequently 
cause damage to the electrodes. An upper limit of 5 A/cm 2 in current density is adopted in this 
design study. 

Since no water-cooled mechanism is planned in the designed MHD generator, the channel wall 
has to serve as heat sink. To prevent burning down of the channel, heat flux through the channel 
has to be limited. For an operation time of 10 to 20 seconds using aluminum oxide as the 
insulator in the MHD channel, heat flux is restricted to below 300 watts/cm”. 


3.3 Generator Performance 

The design calculations are performed by using an ERC’s in-house MHD generator computer 
code, GEN- ID. This computer code is based on a one dimensional generator flow model which 
approximates the influence of the boundary layer in the plasmadynamic calculation. However, 
this computer model does include a direct accounting for the boundary layer in its modeling of 
the generator electrical performance. The local voltage drop across the boundary layer is 
determined by assigning the normalized voltage drop factor, A , defined as V d /(u B h), where V d 
is the effective voltage drop and h is the channel height. For small generator, A is typically 0.3, 14 
this will reduce the generator efficiency by 30% and power density by a factor of 2. 

In GEN- ID, the real gas equation of state is provided by a prior combustion calculations utilizing 
the NASA SP-273 thermo-equilibrium computer program. 12 This program computes the thermo- 
equilibrium process to provide a real gas Mollier chart which is utilized in GEN- ID to determine 
state points through either user defined curve fits or chart interpolation. All the thermodynamic 
variables are fit as functions of pressure and temperature. The combustion calculations include 
solution for the plasma electrical transport properties. 13 Both electrical conductivity and mobility 
are also represented as functions of the local thermodynamic states (pressure and temperature). 
Using GEN- ID, the flow can also be assumed to act like a chemically and vibrationally frozen 
flow. 

Using GEN- ID, parametric studies of generator performance are performed as a significant part 
of generator design. These studies considered parametric variations of inlet Mach number and 
channel geometry while other parameters treated as fixed value variables are listed in Table 2. 
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Mass Flow Rate (kg/s): 


Combustor Pressure 300 psi 

0.228 (potassium) 
0.235 (cesium) 

Combustor Pressure 1000 psi 

0.75 (potassium) 
0.77 (cesium) 

Seeding Percentage (wt.%) 

3.74 (potassium) 
9.66 (cesium) 

Magnetic Field Strength (Tesla) 

1.5 

Effective MHD Channel Length (inch) 

20 

Electrode Pitch (inch) 

1 

Diffuser Exit Pressure (atm) 

1 

Channel Wall Temperature (K) 

600 

Wall Roughness (mm) 

0.03 

Prandtl Number 

0.7 


Table 2. Fixed Value Parameters Applied in Faraday Generator Design 


For each electrode pair, loading factor K, defined as the ratio of external load resistance to the 
sum of external and plasma resistance, is initially set to its optimum value for producing 
maximum electrical power, i.e. 0.5. If the constraint of axial electrical field or that of Faraday 
current density is violated, loading factor will be increased to satisfy the constrain. 

Chemical non-equilibrium effects become significant in high speed nozzle flows. GEN- ID does 
not have the capability of treating these non-equilibrium effects. Instead, for each design case, 
calculations are performed with both chemically equilibrium and frozen flow assumptions and 
their results are compared. Frozen flow calculations are activated downstream to the point 
where local Mach number reaches 2.4 in the nozzle. 

A summary of the findings obtained from the parametric studies are: 

1. An MHD generator with constant velocity was thought to be a good 
configuration to resume the rocket thrust, it turns out to be too much expansion 
in the MHD channel (expansion ratio >10) after certain trials. Constant Mach 
number channel is a better starting point. 

2. The MHD interaction is a strong function of combustor pressure due to plasma 
conductivity level. Cases with combustor pressure of 1000 psi are unable to 
satisfy both constraints of magnet size and heat flux simultaneously. 

3. The generator inlet Mach number has to be high enough to reduce the flow 
temperature to limit heat loss rate. 
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Starting from a constant Mach number channel design. Figure 7 shows the distributions of 
important flow and electrical properties, as well as heat transfer rate in such an MHD channel 
with potassium seed. The channel inlet has a area expansion ratio of 8 with respect to the nozzle 
throat and the combustor pressure is 300 psi. The results show that too much expansion in the 
MHD channel (expansion ratio = 6.5), the pressure drops very fast, as a result, the Hall 
parameter exceeds the design criteria in the downstream portion of channel. This large Hall 
parameter causes strong axial electric field and the loading parameter is adjusted from its 
optimum value to limit the local electric field, which reflects in the plateau in the Hall E-field 
distribution. Also, relative high plasma temperature at the inlet induces a locally high heat flux, 
breaks the design constraint too. Even though these conclusions are obtained under the 
assumption of thermo -chemical equilibrium, the results found by using a frozen flow assumption 
also show the same trend and are thus omitted from listing herein. 

Next, efforts focus on removing those violations in design constrains. In practice, the channel 
contour is linear for simplicity of fabrication and side walls should be designed with a constant 
slope (i.e., constant wall angle). Smaller wall angles are chosen to limit flow expansion within 
the MHD channel. A higher inlet Mach number is selected to locally decrease the flow 
temperature and the associated heat flux value. After all the designing constraints are followed, 
adjustments in the aspect ratio of channel cross section are carried out to provide more 
allowance in the magnet pole gap for the channel wall thickness. Thicker walls provide better 
heat sink effects. During the parametric study, the generator performance and the resultant 
thrust of the rocket based MHD system are of primary concerns. That is to say, the goal of the 
parametric study is to search for an optimum group of design parameters for an MHD generator 
based on an existing rocket motor and magnet. This optimum design parameter combination 
satisfies all the design constraints and, at the same time, produces electric power as much as 
possible and minimizes the influences on the thrust of the rocket system. 

After several adjustments, the designed channel geometry is obtained and included in Figure 8. 
The constant width MHD channel has a square cross section at inlet (1.67 x 1.67 inches) and an 
expansion half angle of 4 degree. At inlet, the area expansion ratio is 9.82 with respect to the 
nozzle throat. 

Figures 9-12 and Table 3 summarize the predicted generator design performance for a potassium 
seeded MHD generator with a 300 psi combustion pressure. Following an expansion flow, 
temperature drops and charged particles (ions and electrons) of a plasma recombine. For an 
expanding plasma, frozen effects prevent recombination and the associated energy release. As a 
result, a frozen plasma flow under expansion reserve higher conductivity but lower plasma 
temperature than an equilibrium expansion flow. This frozen effect can easily be observed when 
comparing data shown in Figure 9 to that given in Figure 10. Stronger MHD interaction, higher 
electrical power, and lower heat loss are induced by frozen effects. 
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igure 7 Calculated Flow and Electrical Parameter Distributions for the Constant Mach Number Faraday Generator. 
RP1 + 02 + KOH, K(wt.%) = 3.74, P comb . = 300 psi, B = 1.5 T, Mass Flow = 0.228 Kg/s, A MH D.mieM,hroai = 8., 
Channel Length = 20 inches. (Equilibrium Flow Assumption) 



1.67“ 



Figure 8 Geometry of the Designed MHD Channel 


In real life design, it is reasonable to assume that the predicted MHD performance of a designed 
generator falls between the gap calculated separately by equilibrium and frozen assumptions. 
Average values of results obtained from these two extremes, equilibrium and frozen, are deemed 
as acceptable design values for the generator geometry proposed in the present design study. Use 
of these average values is supported by comparisons of different MHD parameters given in 
Figures 1 1 and 12, that is, the gap between two extremes of each parameter listed is not wide. 

Therefore, the predicted power output for the proposed MHD generator is 54.2 KW, with 
potassium seed and 300 psi combustor pressure. The power extraction corresponds to an 
enthalpy extraction ratio of 2.36%. The thrust for this rocket based MHD system is estimated to 
be 90 Ibf, assuming ideal expansion to atmosphere after the generator exit. 

Once cesium hydroxide is used as seed component in the proposed rocket based MHD system, 
the theoretically predicted performance are summarized and included in Table 4 and Figures 13- 
lb. Comparison between the results given in Figure 14 to that in Figure 13, relative to the 
equilibrium calculations, stronger MHD interaction of the frozen flow assumption only available 
in the upstream portion of the MHD channel. Too fast flow velocity drop in the upstream portion 
results in relatively lower MHD interaction in he downstream part and the total enthalpy 
extraction ratio of frozen flow assumption is lower than that for equilibrium flow assumption. 

Obviously, due to higher conductivity of cesium plasma, the MHD interaction in the channel 
becomes stronger and the generator performance improves, relative to that obtained with 
potassium seed. The estimated power produced by the same rocket based MHD system using 
cesium is 92 KW, equivalent to an enthalpy extraction ratio of 4.16%, and the thrust becomes 
72.6 lbf when expands to atmosphere isentropically. 
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Frozen 


Equilibrium 


Total Mass Flow Rate 0.228 

Thermal Input 2.24 

Entrance Total Pressure 21 .2 

Exit Total Pressure 6.69 

Entrance Total Temperature 3713 

Exit Total Temperature 3419 

Generator Length 20 

Generator Mean Flow Area 0.00333 

Generator Volume 0.0017 

Internal Surface Area 0.123 

Average Heat Flux -154 

Mean Static Temperature 2448 

Mean Velocity 2639 

Mean Mach Number 3.02 

Mean Conductivity 13.8 

Mean Hall Parameter 2.62 

MHD Power 43.7 

Specific Power Output 0.192 

Interaction Parameter (Ip)* 0.723 

Power Density 25.8 

Enthalpy Extraction 1 .95 

Generator Efficiency 19.9 

Electrical Efficiency 36.5 

Open Circuit E-field 2.77 

Short Circuit Current 20.5 

Internal Impedance 68.5 

Joule Heat Dissipation 42.0 

MHD Push Power -120 


0.228 kg/s 

2.35 MW 

19.8 atm 

4.61 atm 

3623 K 

3290 K 

20 inches 

... 0.00333 m 2 

0.0017 m 3 

0.123 m 2 

... -148 W/cm 2 

2370 K 

2571 m/s 

2.92 

.. 21 .4 Mho/m 

2.74 

64.7 W 

0.283 MW/kg 

1.06 

..38.2 MW/m 3 

2.76 % 

27.5 % 

36.4 % 

2.70 kV/m 

31.0 A 

....44.2 Ohms 

62.2KW 

-178 KW 


ri J X B 

* MHD Interaction Parameter Based Upon Pressure, / = I dx . 

1 J 0 p 

Table 3. Predicted MHD Generator Performance (Potassium Seed, P c ,mb. = 300 psi) 
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M* M/s Atm °K 

0.015 3.1 r ^2800 0.5 r ~i2600 
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Figure 10 Calculated Flow and Electrical Parameter Distributions for the Designed 4° Half Wall Angle Faraday Generator. 

RP1 + 02 + KOH, K(wt.%) = 3.74, P com b- = 300 psi, B = 1.5 T, Mass Flow = 0.228 Kg/s, A MHDil ,i e /A lhroal = 9.82, 
Channel Length = 20 inches. (Frozen Flow Assumption) 



Faraday Current Faraday Voltage 

200 r 
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Figure 1 1 Loading Conditions for the Designed 4° Half Wall Angle Faraday Generator Using Potassium Seed. 

(Solid Lines: Equilibrium Flow, Dashed Lines: Frozen Flow) 
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Frozen 


Equilibrium 


Total Mass Flow Rate 0.235 

Thermal Input 2.14 

Entrance Total Pressure 21 .2 

Exit Total Pressure 3.23 

Entrance Total Temperature 3731 

Exit Total Temperature 3285 

Generator Length 20 

Generator Mean Flow Area 0.00333 

Generator Volume 0.001 7 

Internal Surface Area 0.123 

Average Heat Flux -128 

Mean Static Temperature 2506 

Mean Velocity 2447 

Mean Mach Number 2.85 

Mean Conductivity 32.5 

Mean Hall Parameter 2.39 

MHD Power 89.7 

Specific Power Output 0.382 

Interaction Parameter (Ip)* 1.53 

Power Density 53.0 

Enthalpy Extraction 4.19 

Generator Efficiency 37.5 

Electrical Efficiency 35.0 

Open Circuit E-field 2.57 

Short Circuit Current 44.7 

Internal Impedance 29.2 

Joule Heat Dissipation 90.2 

MHD Push Power -0.256 


.... 0.235 kg/s 

2.28 MW 

1 9.8 atm 

2.83 atm 

3641 K 

3252 K 

20 inches 

...0.00333 m 2 

0.0017 m 3 

0.123 m 2 

.. -147 W/cm 2 

2432 K 

2402 m/s 

2.78 

..35.9 Mho/m 

2.58 

94.0 W 

0.400 MW/kg 

1.64 

. 55.5 MW/m 3 

4.13% 

35.5 % 

35.0 % 

2.52 kV/m 

48.7 A 

... 26.3 Ohms 

97.1KW 

-269 KW 


* MHD Interaction Parameter Based Upon Pressure, 


',=1 


i J xB 


dx . 


Table 4. Predicted MHD Generator Performance (Cesium Seed, P com b. = 300 psi) 
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Figure 13 Calculated Flow and Electrical Parameter Distributions for the Designed 4° Half Wall Angle Faraday Generator. 

RP1 + 02 + CsOH, Cs(wt.%) = 9.66, P comb . = 300 psi, B = 1.5 T, Mass Flow = 0.235 Kg/s, A MHD , inle /A lhroal = 9.82, 
Channel Length = 20 inches (Equilibrium Flow Assumption). 
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Figure 15 Loading Conditions for the Designed 4° Half Wall Angle Faraday Generator Using Cesium Seed. 

(Solid Lines: Equilibrium Flow, Dashed Lines: Frozen Flow) 



550 «- 100 
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4. Conceptual MHD Channel Design 


Because of the limited scope of the proposed testing program and because it is necessary to limit 
the development costs of the MHD channel, we based our conceptual MHD channel design on a 
heat sink configuration. Maximum operating times are not planned to exceed 10 to 15 sec. The 
detailed design strategy is based on practical experience with similar heat-sink configurations. 

A schematic of the proposed electrode wall assembly for the Faraday channel is given in Fig. 17. 
This concept utilizes a phenolic back wall for insetting the electrodes and inter-electrode 
insulators. Phenolic is a plastic material which displays good electrical insulation characteristics 
and reasonably good high temperature resilience. The copper electrodes are attached to the back 
wall using standard fasteners. In addition, the electrodes are machined to accept thin aluminum 
oxide insulators. The recessed cuts on the electrodes holds the insulators in place, and the 
grooves are grouted with aluminum oxide filler. 

The insulating sidewall also uses a phenolic back wall for insetting thermal and electrical 
insulating material. We propose to use aluminum oxide blocks or aluminum oxide casting 
material as the insulating material. The conceptual design for the sidewalls is shown in Fig. 18. 
This method has been proven to provide adequate thermal protection for short length tests. 

The cross-section of the heat sink channel design, indicating the box wall construction method, 
is shown in Fig. 19. Note that the channel is based on a rectangular cross-section having an 
aspect ratio governed by the warm bore air gap of the magnet. 

A schematic of the small-scale rocket based MHD system is shown in Fig. 20. Note that the 
modified rocket nozzle uses an adapter to change from a circular cross-section to a rectangular 
cross-section. 

A schematic of the testing configuration for the small scale rocket based MHD system is shown 
in Fig. 21. This diagram indicates the propellant feed systems, auxiliary power systems, and 
measurement and control systems needed to test the rocket based MHD concepts. 
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Electrode Wall Assembly 


Phenolic Back Wall 



Copper Electrode 



Figure 17 Conceptual Design of Electrode Wall Assembly for Faraday Channel 
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INSULATING WALL 
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HEAT SINK CHANNEL ASSEMBLY 
Phenolic - Box Wall Construction 

Figure 19 Conceptual Design of Heat Sink Faraday Channel Assembly 
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Magnet 



MHD Channel 


Figure 20 Schematic of Small-Scale Rocket Based MHD System 


LIQUID ROCKET MHD GENERATOR TEST APPARATUS 



Figure 21 Schematic of Testing Configuration for Small-Scale Rocket Based MHD System 
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5. Rocket Based MHD Accelerator 


Detailed descriptions of the potentials of MHD accelerator technologies applied in hypersonic 
flight and sea water propulsion were given in References 15, 16, and 17, respectively. Following 
is a summary of the operation of an MHD accelerator and an estimate of performance potential of 
an MHD accelerator using operation conditions similar to that coupled to a rocket motor. 


The MHD accelerator is composite system which consists of an MHD channel, a magnet, and a 
power supply as its major features. Operation of the accelerator consists of establishing the flow 
of a conducting fluid through the MHD channel, which is situated within a magnetic field 
directed normal to the flow. By application of an external power supply over the channel, an 
electrical current is forced across the conducting fluid. The action of Lorentz force on the fluid 
will then accelerate it, this action is manifested in the form of an increase in fluid total pressure. 
With the MHD accelerator channel designed for proper flow expansion through its length, the 
total pressure increase is in the form of velocity, or alternatively, flow kinetic energy. Joulean 
heat dissipation increases the total enthalpy of the flow. That is, within the MHD accelerator, 
simultaneous increase in flow total enthalpy and total pressure occur through the effects of 
Lorentz force as well as Joule heating. It is this fact which is often quoted when referring to the 
benefits of MHD flow acceleration over other heating technologies. Direct heating a channel 
flow increases the total enthalpy of the flow but, on the other hand, reduces the total pressure. 
The latter represents a loss to the effectiveness of the acceleration process. 


The ideal Faraday accelerator is conceptually one with infinitely electrode segmentation such that 
no Hall current is present in the plasma. As a result, the Lorentz force is aligned with the flow 
direction and no side forces resulted. This is the type of MHD accelerator examined in this study 
since creation of uncontrollable side forces in the flow is generally unfavorable during the 
maneuver a rocket based system. Formulation of electrodynamics in the channel of a Faraday 
accelerator can be easily derived based on Ohm’s law as 


axial current 

j, = o 

(5) 

Faraday current 

j y - c u B (K -\) 

(6) 

axial field 

E=onuB(K-\) 

( 7 ) 

Faraday field 

E y = K uB 

(S) 


where the fields and currents are taken as integrated averages over the cross section, and K is the 
normal power supply load factor defined as 




K = 


applied 


V 


>pen~ circuit 


uB’ 


K> 1 


(9) 
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For a Faraday accelerator, the electrical efficiency can be written in terms of fundamental plasma 
variables by substitution in Equation (9). That is. 


n 



( 10 ) 


The expression (10) for the electrical efficiency shows that the most efficient Faraday accelerator 
is one operated with low load factor. 

The Lorentz force and power supplied terms are then written as 

u(JxB) = au 2 B 2 {K- 1) (11) 

J E = y v £,. = <7 u 2 B 1 K{K- 1) (12) 


The above expressions emphasize the need to maximize both the magnetic field and conductivity 
to enhance MHD interaction. The general performance of an ideal Faraday accelerator with a 
conductivity of 30 Siemens/m is depicted by the plots of Figure 22. These graphs illustrate the 
thrust current, electrical field, input power required, and push work produced for lines of constant 
induced electrical field («x5). For any given velocity and magnetic field combination within 
the range shown, the power density requirement can be estimated. 

For example, for a rocket based MHD accelerator operating at a load factor of 1.1 {T] = 90%) 
with a global plasma conductivity of 30 Siemens/m, velocity of 2500 m/s and a magnetic field of 
2 Tesla ( uB = 5 kV/m). At this condition, the ideal performance (no flow field losses) of the 

'J 

accelerator found from Figure 22 is current density 1.7 Amps/cm“, applied field needed is 5.6 
kV/m, 93 MW of electrical power for every is required, and 83 MW of MHD work per cubic 
meter of volume is produced to accelerate the plasma flow. Increasing magnetic field for the 
same operating point (moving to higher u x B lines) increases efficiency. At any given uxB 
value, higher power supplied or current density, leads to lower efficiency. As a result, to 
maximize efficiency over a specified speed range, low current and high magnet field are 
required. 

The effects of increasing conductivity are shown in Figure 23 where increasing conductivity, for 
the same ux B value, decreases power requirement, promotes push power and leads to more 
efficient operation. 

In consideration of a constraint on current density of 5 Amps/cm" as proposed in the generator 
design, the efficiency obtained from Figure 22 has to be no less than 75, 86, and 90% when uB 
has values of 5, 10, and 15 KV/m, respectively. Under such circumstances, the upper limits of 
the push power density are relatively 250, 431, and 750 MW/m 3 for uB of 5, 10, and 15 KV/m, 
while the corresponding limits on the required power are 333, 501, and 833, respectively. 
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Figure 22 Performance Diagram for Ideal Faraday MHD Accelerator with a Conductivity of 30 Siemens/m 

(Solid Lines: uB = 15 kV/m, Dashed Lines: uB = 10 kV/m, Dash-Dotted Lines: uB = 5 kV/m) 



Input Electric Power MHD Push Power 




Figure 23 Comparison of Performance for Ideal Faraday Accelerator with Different 
Conductivities 


6. Conclusions and Recommendations 


A Faraday type of MHD generator driven by a sub-scale LOx/RP rocket motor is conceptually 
designed based on an existing magnet and by assuming realistic limits on (1) the axial electric 
field, (2) the Hall parameter, and (3) current density, and (4) heat flux (given the criteria of heat 
sink operation). A parametric study is performed using GEN- ID computer code to find the 
optimum geometry of the MHD channel which not only satisfies all the specified constraints but 
also maximizes the electrical power. This study leads to a constant width MHD channel with a 
diverging half wall angle of 4 degree based on a combustion pressure of 300 psi. With a 
magnetic field of 1.5 Tesla, the electrical power output from this generator channel is estimated 
to be 54.2 KW with potassium seed and 92 KW with cesium seed. The parametric study shows 
that cases with 1000 psi combustion pressure can not simultaneously fit in the existing magnet 
structure and satisfy the designated heat flux constraint. 

A conceptual design of the Faraday MHD channel and a testing configuration for the entire sub- 
scale rocket based MHD system are proposed. This Faraday Channel utilizes a phenolic back 
wall for inserting the electrodes and inter-electrode insulators. Copper electrode and aluminum 
oxide insulator are suggested for this channel. The thickness of the insulator wall is determined 
for a maximum operating time of 10 to 15 seconds. 

The present design study is totally based on the ECF magnet and assumes that it can be 
refurbished to produce a magnet field strength of 1.5 Tesla. To fit into the narrow gap of this 
existing magnet, results of the parametric study show that the MHD channel must have a 
rectangular cross section to provide the required mass flow rate and to leave enough allowance 
for the wall thickness required for heat sink operation. A transition nozzle from circular to 
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rectangular cross section is thus needed to connect the rocket combustor and the MHD channel. 
Therefore, in addition to using this low cost magnet, it might be worthy, at the same time, 
searching for other candidate magnets which can provide higher magnetic strength, enough 
space for the MHD channel expansion, and consideration of employing a channel with circular 
cross section. Matching the circular rocket combustor with a circular MHD channel will result 
in a channel design much easier to be fabricated. 

In this study, chemical non-equilibrium effects in high speed flow are accounted for by 
including certain degree of frozen flow effects to a chemically equilibrium model. To more 
precisely predict the influence of the non-equilibrium effects, modeling work including kinetics 
of different reactions should be developed. Also, multiple dimensional model is needed for 
better estimate of the boundary layer effects and the side forces effects in the MHD channels 
other than Faraday types. 

It is strongly recommended to continue the rocket based MHD generator research program to 
testing phase in order to demonstrate power production from the designed generator, to set up a 
data base for verification of analytical modeling work, and to accumulate experience in 
designing and operating an MHD generator of this kind. Once these goals are reached, large 
scale rocket based MHD generators can be later designed to provide great performance with 
confidence. 

Performance of a rocket based MHD accelerator is estimated without consideration of flow field 
losses. In this simple approach, current density, push power density, applied electrical field, and 
input power density can be found when conductivity, loading factor (or efficiency), and induced 
electrical field u x B are assigned. The results show that the lower the MHD interaction, the 
higher the efficiency in an accelerator. Higher conductivity decreases power requirement, 
promotes push power, and leads to more efficient operation. It is also found that with a current 
density constraint of 5 Amps/cm 2 and a conductivity of 30 Siemens/m, the push power density 
can be 250, 431, and 750 MW/m 3 when the induced voltage uB have values of 5, 10, and 15 
KV/m, respectively. 

More thorough analytical modeling of this rocket based MHD accelerator system taking into 
consideration of flow field losses is deemed absolutely necessary to better estimate the 
performance potential of this type of MHD accelerator under other practical constraints. 
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APPENDIX A 


REVIEW OF RUSSIAN WORKS 
ROCKET BASED MHD TECHNOLOGY 

Prof. Valentin Bityurin 


As part of this Task 1, ERC consigned its Russian consultant. Professor Valentin Bityurin 
to provide a synopsis on past and current Russian technology related to rocket based MHD 
concepts. 

This Appendix presents the results of that background study. 
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MHD GENERATORS BLTLT-IN IN ROCKET ENGINES. 
STATUS AND PROSPECTIVE 

(Review of Russian Works) 


3 8 
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Abstract 


In this report scientific information on R&D works available from open literature 
of Russian specialists on MHD built-in the of a rocket nozzle engine are present. It 
seems that the pick of activity in this field is passed 5-10 yeas ago for very clear reasons. 
It should be noted also that experimental investigations made are mostly occurred out 
of publications. 
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INTRODUCTION 


Much attention has recently been given to the concept of all-electric aircraft [1.2]. 
To improve the tactical and technical characteristics of such vehicles, it is proposed to 
equip them with novel electrical and power facilities. Atmospheric vehicles with gas- 
turbine propulsion units utilising relatively low-temperature flows of a working medium 
(T<2000K) are preferably provided with high-power electric sources such as 
electromechanical (in particular, synchronous) generators driven by the main 
propulsion units of the aircraft [3.4]. From the other hand it appears promising to equip 
aerospace vehicles with high-temperature (T>3000K.) rocket engines with MHD 
generators (MHDG) used as high-power airborne power sources [5-7], The use of 
MHDG in the power facilities will make it possible to produce new types of electrical jet 
engines for aircraft, replace the conventional turbopump units in liquid-propellant 
rocket engines by advanced fuel supply systems, and accomplish the MHD control of 
the thrust vector at aerospace vehicles [8-12]. An MHD generator installed in this nozzle 
of a scramjet is considered as a one of the main component of an deep optimized thrust 
system of hypersonic transatmospheric aircraft [12]. 

Research and development works in field of interplanet voyages and. in particular, 
a manned mission to the Mars considers a nuclear power production/thruster facility 
based on a nuclear gas-phase reactor with magnetic stabilization of active material as 
one of the most prospective solution [13.15], For such an on-board power generating 
facility the MHD generator located somewhere in nozzle of the nuclear rocket engine 
could be a rather effective electricity generating unit [15,16]. 

In the most papers devoted to airborne high-power MHD generators (self- 
contained or built into the nozzle of the propulsion unit). MHD channels are discussed 
that are based on a two-pole scheme. The processes of energy conversion in these MHD 
generators are usually accompanied by disturbances of flow in the entire channel 
volume. Being an auxiliary power generating most built-in MHD generator should be a 
low disturbing device that does not decrease significantly the overall performance of the 
power generating system. For this reason the MHD generator working volume should 
presumably occupy only small portion of the whole working volume of the power 
facility. One of the approach helping to solve this problem is an utilization a multipolar 
magnetic system instead of the traditional two-pole system. Depending on the number 
of pairs of poles, the proposed facility makes it possible to accomplish MHD 
conversion both in the nozzle peripheral zone and in the central part of the flow. 
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SCHEMES OF MULTIPOLAR MHDG 


In the multipolar MHD generators under consideration [6.7], the electrodes may 
occupy a part of the surface of thin cylindrical rods (Figs. I a and lb) or be installed on 
U-shaped projections (Fig.lc) and on the side walls of the channel (Fig. Id) [6.7]. The 
magnetic Field B is produced by a multipolar magnet system (MS), which consists of 2 P 
current-carrying linear sections located above the electrodes. With P> 2. such a magnet 
system has low scattering fields as compared to the traditional magnet system, and this 
is important for the airborne multipolar MHDG. The current-conducting gas moving 
with velocity v and interacting with the magnetic Field causes electric current in the 
load. The external wiring of the electrodes of a multipolar DC MHD generator and the 
definition of basic geometric parameters of MHD channel having modules of different 
are shown in Fig. 1 . 

In the other versions electrodes are set on the U-shape projections (Fig. l.d) or on 
the cilindrical channel walls (Fig.l.e). The external magnetic Field B is formed bv' a 
multipolar magnetic system containing 2P linear sections with the transverse dimensions 
/■\( <</u. Outside the channel the sections are located on radius r\\>ri> directly above the 
electrodes. Fig.l.f shows the external electrode commutation for DC and one-phase AC 
generators. It is aiso possible to get a three-phase AC from MHDG if the magnetic 
system is excited by a three-phase AC and another versions of the electrodes 
commutation is used. 



THEORETICAL MODELS 


Further it assumed anways that the magnetic fields induced by the electric currents 
in MHD channel can by neglected. The period of AC is much more than the period of 
the gas flight through the channel. The length of MHDG L is considered to be much 
large than the distance between the electrodes. The end effects in channel are rather 
small. The supersonic and subsonic flow regimes at the finite values of Hall parameter 
([3o^0) and high hydrodynamic Reynolds numbers (Re>>l) are studied. This allows to 
describe three-dimensional processes in the channel usual MHD approximations for a 
stationary nonviscous ideal gas flow [17,19], (Additionally the influence of the thin 
boundary layers is taken into account when the electrodynamic problem is further 
discussed.). The dimensionless MHD equations are the following: 


div(p u.)=0: (1) 

p.(u.grad)u.=-( l/(yM) : )grad p.+S/XJ-B: (2) 

p. T(u grad)j.=SM- 7 ( 7 - 1 )J. : /(/- 0 .): (3) 

J ,-r(JW|B|)J xB=a.(uxB.+E.): (4) 

div J =0: E=-gradC: (5). (6) 

div B.=0: rot B.=0 (7). (8) 


p.-p.T.: s=ln(/?./p.): a=G.(p.. Ty. (3o=f3eo(/L. T.)|B|; u=u./uo: 

T = T.iTu p=pjp„: p.=pi/po: cr.=at/'Oo: E.=EJ(uqBo): 

J — J./( o tjinBo ) : B—B./B o. 

Here the inferior index «</» corresponds to dimensional values: «0» — to base 
values at the inlet of MHDG; «*» — to dimensional values: // [u^.u^u.) is a gas velocity: 
p. p , T and v are a density, pressure, temperature and gas entropy respectively: 

Gp =G () /(l+ ) is an effective plasma conductivity: J. E and U are an electric current 

density, electric field and electric potential respectively: M = u 9 /^7 R 71) and 

S = Gj3 0 5fi/(po«o ) ars nutT >ber a parameter of MHD interaction: y. R are an 

adiabatic coefficient and gas constant: X=L/r„ is a relative channel length: 


B t) = '/fjcPm ) is a characteristic value of the magnetic flux density 

(it(i=4^xl 0 : ); Iwu=pL is the total current through p sections of MS. 

At a weak MHD interaction (S<<1) the problem is solved by means of 
perturbation techniques, when .dependent variables, such as ?/.. T and etc., are 
presented as the asymptotic expansion series: 

//.=//+S-w ? +ar «»+...: p=p+S-p^&ip(i-r.... (9) 


Here S<<! and £=/;//»<<l are small parameters: the values without index belong 
to zero-order approximation: lower «.y» and «e» indexes correspond to the first 
approximation of perturbanced parameters according to .v and £ respectively. At ,v<<1 
the unperturbanced gasdynamic parameters in MHDG are known. With this 
assumption the values B and J can be found from the independent solution of two 
elliptic problems Eq (4)-(6) and Eq (7)-(8). 

For long cylindrical MHD channels (k=Llr {) > 4 to 5. see Fig. 1) the distribution of 
B and J away from the end zones is found from the independent solution of two elliptic 
problems of the form [19] 


AA r = 0 
B = rot( .4. . i - ) 


A r(D 'E-A, o ln[a/(l-p-)jV^ = 0 

J = a[v,^-i ; 3 V r(0 4»]/(l-(3 ; ) 


( 10 ) 


where A(O.O.H-) and v F=vA -U are the vector potential of the magnetic field and the 
electric current density potential, respectively: U is the electric field potential in the 
MHD channel: a and v are the electrical conductivity and velocity of gas. respectively. 

After that one can find perturbanced fields of the corresponding gasdynamic 
values (n,. /?... p ) from the corresponding solution for M>1 of the inhomogeneous 
hyperbolic problem of the form [19]. 


v 5 =-grad nK x + J (JxB)«fc 

o 


(ID 


where % is the velocity disturbance potential v y r.): and /(J.B.r) is the known 

function of J.B.r. whose structure is determined by the design of the electrode modules 
(Fig. 1). 
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At a high MHD interaction (,v>0.5) the characteristics of MHDG with the variable 
cross-section are determined on the base of one-dimensional MHD equation [17.19], 

The analytical solutions for the calculation of the parameters A.. T. % of 
multipolar MHDG at , v <<1 are given in [19], Figures 2a and 2b show the calculation 
results of the typical distribution of the function of the magnetic flow A. and current 
density potential T' in cylindrical MHDG with U-shaped electrodes. The function 
,-( : =const is given for the case of P - 3 and P=5. The isolines H'jy. tpA,(l - k/k nl )] = const 
are constructed for P= 3. the load parameter A'=0.5 and the typical values of the short- 
circuit current drop factor k\ and no-load voltage factor A'„i (the factors /m and k a i allow 
for the effect of a and r on 'F in the electrode boundary layer). Figure 2c also shows the 
distributions of the dimensionless velocity disturbances y ? (v jr . y S v . r._.) in an MHDG with 
\=L : Jr<\= 3 for different Hall numbers ([3=0 and [3= 1 .72). The calculations of A.. T. and 
v . were performed for the case of small cross-sectional size of the linear portions of the 
magnetic system. From Fig. 2. it is clear that the distributions of ,-f. T. and r r at P > 2 are 
nonuniform in the channel cross section. The process of MHD energy conversion is 
localised in the peripheral zone. At (3^0. the axial symmetry of the disturbance fields of 
the velocity r . is broken, and cross eddy currents are formed in the MHD channel 
outside of the flow. This effect results in deviation of disturbed flow to the cathode 
electrodes. 

It is well known that ideal single-pole VIHDG (P- 1) with electrodes arranged on 
the channel side walls (as shown in Fig. Id) feature higher energy characteristics at a 
constant level of the magnetic flux in the interelectrode zone. However, their technical 
realisation in airborne power facilities has a number of disadvantages. These include the 
contraction of currents on the electrode edges, high scattering fields of the magnet 
system, high consumption of ionising seed, disturbances in the entire nozzle volume 
[19]. As mentioned above, multipolar MHDG with U-shaped electrodes are free from 
these disadvantages and provide diffusion current flow in the channel, low scattering 
fields and a low flow rate of the ionising seed in the zone of effective MHD-interaction 
[ 6 ], 

The output parameters of multipolar MHDG with s< 1 in MHD channels that are 
barely expanding along their length can be estimated using quasi-one-dimensional 
dimensionless MHD-equations. averaged over the channel cross section, of the form 
[19] 
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( 12 ) 


p'- : n.i i — 


710 




'o 


pr 


dv - 1 d P 1 / T „\ j c ; C ;> L : 

- ; + .V — J X B). -r 21 -r pr. 

dz '/VI - dz A ; D D„ , ' 


(13) 


d 


pv. -f [ r - V : M : 1 v : |= 5 M"(y - 1)(J B )±+2C t ^(T,- T h ) 

a* l _ ! A ” U f 


. 1 


pv. 


(14) 


c/5 6 c/p 5 c/v. C • 

+ - (2- H) - = — . 

dz p dz v. dz 2 

V - 1 , 

P= p T . T, = T -r 0.9 M ' v: . (15) 


r 0 s - 


Bjj Lu 
Po v o 



Here. C, and C. are the coefficients of friction in the channel and of pressure loss 
on the electrode modules [20.21]: £>.=4/717 and D p - 4 FLi {IPzirp ) are the hydraulic 
diameters: a is the factor of shading of flow by the electrode modules: o* and 5** 
denote the displacement and momentum thickness, respectively: y is the adiabatic 
coefficient: k=£/r.i is the relative length of the MHD channel: p. p. T. and v : . are. 
respectively, the dimensionless density, pressure, temperature, and axial velocity of the 
gas flow: s and M are. respectively, the MHD-interaction parameter and the Mach 
number, determined by the gas parameters at the inlet: H=5*/ o** is the form parameter 
of the boundary layer: 7], is the channel wall temperature. 

Within this approach, the friction coefficient C. is determined on the basis of the 
limiting laws of friction and heat transfer [20]. The cross section-averaged specific 
electromagnetic <JxB>. and electric power <J E> in (12)-(15) are found from the 
solution of the two-dimensional electrodynamics problems on transverse edge effect in 
cylindrical multipolar MHD channels [19]. 

Figure 3 shows typical dependences of the values of <JxB> : . <J E> short-circuit 
current /, c and open-circuit voltage Uni on the size r ; of electrode modules of MHD 
channels with L'-shaped (Fig. 3a) and segmented [Fig. 3b. where 5=(<pi>-(f>i)/(po] electrodes. 
The dots in Fig. 3a indicate the results of the experimental study of the characteristics of 
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multipolar MHD channels [22]. Experimental data have been obtained for E*— 1.2.3 in 
plasma MHD-models with the channel length L c i,= 1.75m and radius /•. ) =0;6m with the 
following parameters of argon flow: flow rate G-2g/s. Hall parameter [3-1-2. effective 
conductivity in the core a,~50 to 100S. MHD-interaction parameter 5-0. 1 [22.23]. From 
Fig. 3 it is clear that the experimental data agree with the theoretical dependences, and 
this fact enables one to recommend them for the calculation of concrete versions of 
propulsion units with a multipolar MHD generator. 

For the given geometry of the engine nozzle and ,y<<l. the approximate analytical 
solutions of system ( 1 2)-( 1 5) for the dimensional quantities of current ! itn and power 
.V ro i of MHDG with U-shaped electrodes have the form [19] 


/, en = 2 tuP J ] 


i A 

A, 

/ V M 

1: fx F ; ! 

Ql 

" 

j -b ] 

1 -J 




! u 0 / 0 ( r n T Pr, [ 

k(z) — U o^ni '' 4v( ’ ) — | ! 

* Ip,; >o [ 



( 18 ) 


Here. U is the MHDG voltage: /<(:) is the local load factor: L is the total current 
(magnetising force) of the magnet system section: Un=47ux I ()-' is the magnetic 
permeability. The typical dimensions of the elements in the A1HDG cross section are 
shown in Fig. I . 


Relations ( 1 6)-( 1 8) make it possible to calculate the output characteristics of the 
built-in multipolar MHDG at v<<! and can be used for optimisation of the power 
facility parameters. 
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EXPERIMENTAL INVESTIGATIONS 


Experimental facilities. The experimental investigations of multipole MHDG were 
carried out on the plasma-vacuum equipment [24j. Argon plasma parameters at the 
inlet of the experimental MHDG had the following values: velocity ~50(H-2000m/sec; 
Mach number M~0. 5^-2. 0: gas temperature ~(8. . . 1 2)xl0 3 K: static pressure ~10 3 Pa: mass 
flow rate G-2-og/sec: average Hall parameter (3- 1-*-2: average effective electrical 
conductivity a-(2-^4)x 10 3 s/m: parameter of MHD interaction S-0.5-G.O. Magnetic flux 
density in the channel was changed from 0.1 to OAT. In experiments (see Fig. 4:) argon 
entered the heater (plasmatron) 1 and then came through the nozzle and MHD-channel 
2 to the vacuum cavity 3. Fig. 4 shows a general view of experimental DC multipole 
MHDG with a cylindrical channel with radius /' (J =60mm. 2P electrodes 4 with length 
L k - 175mm are symmetrically located on radius p c =46mm in the peripheral channel 
zone. The electrodes were made as thin graphite cylindrical bars by radius r e =9mm. The 
external magnetic field in MHD channel was formed by multipole magnetic system 
(MS) 5. 

MHD characteristics. Fig. 5 shows the voltage current UfJ) and power .V,(7) 
characteristics of MHD channel at the different pole number 5=1. 2. 3. In the 
experiments we have an electrical current for MS 7,,„=400A (5-0. 3T). mass fow 
Gi=2g/sec. gas velocity ;/.=400-500m/sec. The dependencies Ufl) for multipole MHDG 
have the falling character. The value of the highest electrical power is obtained in the 
channel at P=2. It can be explained by the fact that at P=2 we have the greatest values 
of the magnetic flux density. Fig. 5 also represents the comparison of dimensionless 
theoretical and the experimental data [6]. It's evident that the experimental data 
qualitatively correspond to theoretical curves. 
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DESIGN AND PARAMETERS OF PROPULSION UNIT WITH MHDG 

Figure 6 illustrates a version of a propulsion unit with a multipolar. MHDG built 
in the nozzle of a liquid-propulsion rocket engine. The facility utilises cryogenic fuel 
components H : ; +0 : The readilvionizable cesium seed is fed directly into the wall zone 
of the nozzle. The relative flow rate of cesium y ( ~G c /G E ~0.3/(P+ 1 ) (here. 0, and G E are. 
respectively, the cesium and fuel flow rates in the propulsion unit) is selected so as to 
provide the required level of ionisation and conductivity in the annular zone A of the 
multipolar MHDG in the corresponding region of effective MHD interaction of the 
flow (A~ro/(P+\)) [19]. The body of U-shaped cooled electrode modules I of the MHDG 
is made of copper and faced with segments 2 of pvrographite on the side of high- 
temperature flow (see Fig. 6). The internal copper-faced surface of the graphite segments 
is soldered to the copper body of the electrodes. The electric current from the electrode 
modules flows through tubular copper current conductors 3 insulated from the body 
and used for the delivery of the cooling agent (Hz:). From the external side of the 
propulsion unit nozzle, the electrodes of like polarity are alternately connected to 
cooled annular current conductors 4. In order to provide electric insulation, the hot 
walls of the nozzle 5 are separated from the electrodes by refractory spacers 6 made of 
concrete with a filler of the MgO type. In the interelectrode zone, the nozzle walls are 
covered with thin insulating modules based on AlzO>. Intensive cooling of the nozzle by 
the oxidizer ( Oz) provides additional protection against electric breakdown thanks to 
the formation of a cold wall boundary layer. The multipolar magnet system consists of 
IP racetrack coils or sections 7. Each section is wound by reinforced wire of ultra-pure 
aluminium .-I9999 having a high electrical conductivity [23]. The sections can be cooled by 
iiquid hydrogen through a system of longitudinal axial channels. The multipolar 
magnet system with power elements is mounted on the external wall of the nozzle in an 
airtight welded housing 8. To remove the power loads of the magnet system arising 
during the MHDG operation, the magnet system housing is additionally secured to the 
propulsion unit carrier frame. The bondage and rigidity of the winding and bandage 
elements of the magnet system were estimated using the methods described in [25]. 

The main parameters of the advanced power facility with a liquid-propulsion 
rocket engine and a built-in multipolar MHDG are given in the table. The properties of 
the products of combustion of the working medium Hzj+Oi ;, with the air excess factor 
(a=0.8. seeded with alkali metal (cesium), and the specific performance of the liquid- 
propellant rocket engine were borrowed from [26.27], The calculation of the output 

characteristics of the MHDG according to (16)-(18) was performed for the pairs of 
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poles of P=3.5. The coefficients of decrease of the short-circuit current k, and the open 
circuit current k ni were taken equal to 0.9 [19.29]. The load factor /<= 0.75. Given in the 
table are the results of the calculation of the parameters of a built-in MHDG for two 
values of the pairs of poles P= 3 and P-5, whose magnet system is characterized by 
weak scattering fields in she far zone. The size of the electrode modules r c =rjn,= 0.05 is 
selected taking into account the fact that the loss of specific thrust of the propulsion 
unit due to the electrodes protruding into the flow is low compared to the losses due to 

MHD-processes in the propulsion unit nozzle (r e « \Jn7c / kGv~ PC., . where C*=l 
[ 1 8]). The dependences of the output parameters for other values of i\, and P are given in 
Fig. 7. 

The analysis of the theoretical data obtained (see the table and Fig. 7) has shown 
that the proposed version of the propulsion unit with a thrust P~X 0.0 ton and a 
multipolar MHDG makes it possible to obtain electrical power A-10MW with low- 
losses of the thrust / A/ > /P-AF’ S p// 5 S p~0.005). With the same magnetizing force (total 
current) in the magnet system, an increase of P at r e P<0.5 causes an increase of the 
useful generator power and a decrease of the size of the annular zone X~r„/P. where the 
disturbance of gas flow is appreciable. The power of a multipolar MHDG may be 
increased at r,P<0.5 w'ith the same composition of the combustion products by- 
increasing the size of the electrode modules. In this case, however, the propulsion unit 
thrust losses will increase due to aerodynamic losses in the nozzle. The power of an 
MHDG utilizing the working medium Hz.+Oz; with alkali seed may be increased by 
adding aluminium powder into the fuel, this increasing the combustion temperature 
and. therefore, the electrical conductivity of gas in the MHD channel [29], 

The proposed scheme with a multipolar MHDG can also be used in power 
facilities utilizing other fuel pairs [e.g.. Oz, + kerosene w-ith a small amount of ionizing 
seeds of an alkali metal (approximately 1-2% of potassium)]. Estimates indicate that the 
effective conductivity of the combustion products of fuels based on Ozj + kerosene 
seeded with alkali metals is of the same order as that for the Ozj+Hzj+C . y fuels [29]. This 
fact allows one to conclude that, in this case, the electrical power of multipolar MHDG 
with the same size of the MHD channel and magnetizing force of the magnet system 
B~\. 5T) will be close in magnitude to the respective theoretical values of .V given in the 
table and in Fig. 7. Note that in propulsion units with MHDG utilizing Oz; + kerosene, 
a superconducting version of the magnet system is preferable. Prior to starting, this 
system is cooled with liquid helium on the ground. The loss of liquid helium caused by 
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the heat flows to the cold zone of the superconducting magnet system on board the 
aircraft will be relatively low (less than I l/h) in the presence of intermediate crvostatic 


screens cooled by liquid oxygen [4]. 

Power facility parameters 

Value 

Table 1 

1. Propulsion unit 

Total flow rate of the working medium G z , kg/'s: 

200 


Flow rate of fuel //^/oxidizer O 2 Gh/G 0 . kg/s 

30/170 


Pressure in the combustion chamber/?.. MPa 

5 


Temperature in the combustion chamber T c . K 

3500 


Expansion ratio with respect to pressure e 

100 


Specific thrust (without MHDG) /%. m/s 

4125 


Engine thrust ignoring. MHDG P.. ton 

82.5 


Thermal power of the propulsion unit V. MW 

2000 


II. MHD channel 

Inlet Mach number M<* 

2.4 


Inlet temperature Tci,. K 

2600 


Inlet pressure P c h, MPa 

0.37 


Inlet conductivity cr C h. S/m 

15.0 


Number of pole pairs P 

3 

5 

Hall parameter p 

0.17 

0.25 

Generator power .Von. MW 

7.0 

10 

Voltage U. kV 

1.6 

1.18 

Electric current I. kA 

4.45 

8.47 

III. Magnet system 

CuiTent density J m . MA/m : 

100 


Magnetizing force MA 

0.75 


Number of sections P 

3 


Induction on the external radius of MHD channel B c 1 ,. T 

1 

1.5 

IV. Weight and size of the power facility 

Critical section radius R cr . m 

0. 1 75 


Average radius of MHD channel R c 1 ,. m 

0.45 


Average radius of magnet system R m , m 

0.52 


Nozzle radius at the outlet R,. m 

0.65 


MHD channel length / c u, m 

0.6 


Magnet system length /,„. m 

0.75 


Total length of the power facility 

(combustion chamber + nozzle) L. m 

2.75 


Total weight of MHDG v , . kg 

200 


Specific weight of MHDG 7 . kg/kW 

0.036 

0.02 

Specific thrust loss of propulsion unit + MHDG: 

P, V mild 2 * P >p 

0.004 

0.006 
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.VIHD CONTROL ESTIMATION OF SUPERSONIC BOUNDARY LAYER 


At the present time MHD control possibility of supersonic boundary layer in the 
problems of internal and external ballistics is widely discussed. It is supposed to change 
pressure distribution and. consequently, balance of aerodynamics forces along surface 
flowing around with use of the decelerating or accelerating VIHD force which able to 
cause a separation or reattachment of supersonic boundary layer [30]. The analysis of 
the turbulent boundary layer separation carried out in our paper [30] has been based on 
the traditional approach widely used for example in paper [31] to treat separation 
phenomena in conventional gasdvnamics. 

A positive pressure gradient in the flow leads to the flow deceleration in the 
inviscous core as well as in the boundary layer. Due to a difference between the core 
and boundary layer velocity (the boundary layer velocity is significantly lower) positive 
pressure gradient results in a much higher deceleration near the wall then in the flow 
core. Thus, the velocity profile becomes less filled and the friction tension 


C f 


9 u e 


at the wall decreases. In order to illustrate the above scheme generalized experimental 
data in terms of the friction coefficient and the velocity profile parameter n (w/w ; =(v/8) l/ " 
as the function of normalized pressure gradient 



are plotted in Figs. 8 and 9. These values are 


0.0128 1 
Re 0 . a ” 1 t 0.22G] ’ 


(19) 


'.266 y l + G | 5/ 9, 
1+0J6G 


where 

0, = j°— (1 — — )d\ — kinematics thickness of momentum loss: 
' f) w u e 


5 1 


ERC-R-97-01 7 


Re 


e. 


p« ,0, 
u 


— Reynolds Number based on momentum loss thickness 9, and 
characteristic boundary layer temperature: 


<5/. <5 — boundary layer and viscous sublayer thickness, 
index «e» corresponds to external border of the turbulent boundary layer. 

At the separation point C,=0 and from the friction law (19) the separation criteria 
can be derived in the form 


Re 0 a " -^-— = 0.058 . (21) 

p«; dx 

In Fig. 10 the comparison of this criteria with experimental data from paper [3] are 
presented in term of value 

8 dp 

P e u] dx 


The flow separation under external counter pressure gradient develops through 
the external pressure penetration into the nozzle along the subsonic portion of the 
boundary layer. The gas flow deceleration caused by such a pressure gradient results in 
a strong growth of the boundary layer. Streamlines are declined apart the nozzle walls 
providing the compression in the inviscouse core flow. Intersection of the compression 
waves creates an oblique shock in inviscous core and the pressure growth till the 
atmospheric pressure that provides the external flow and boundary layer flow 
conjunction. The process of viscous-inviscouse interaction on the external pressure 
penetration length is described approximately by the separation criteria (21) and 
Prandtl-Mayer expression for pressure growth in a week compression wave. 

dp 1 d8 

Pc“r~ V .V/, : -! dx 


Then the dependencies of the critical pressure ratio pip i and the interaction zone 
size /, on flow parameters with empirical coefficients derived from experimental data 
comparison are as follow 


El 

Pi 


'P-5 

4- 0.6841 — J Re, 

Pc ) 


- 0.1 


JM/ 


( M~ - 1 ) 


0.25 


( 22 ) 
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(23) 


/., = 19.10. y ; Re 0 OI25 ( A7, : - 1) -0 - 25 . 

It was found that the above defined values for different cases of supersonic flow 
separation: in front of step obstacle, a wedge, under incident shock wave interacting 
with boundary layer and under off-design nozzle flow are practically the same. It means 
that the process is defined by the selfsimilarity and. consequently, flow parameters 
distribution in the separation point region as well as size of the latter are independent 
on the reason of the separation. This fact will be used to derive the separation criteria 
for MHD interaction case. 

Despite the great practical interest the turbulent layer separation phenomena in 
MHD generator flows are studies insufficient up to now. For example, there are many 
publications on experimental results about MHD separated flow in «PAMIR» type 
MHDG [17.32-34] and only one paper by A. V. Gubarev [35] devoted the more or less 
closed theoretical model of MHD separated flow description. 


MHD Separation Criteria. 


Let us consider an incompressible turbulent boundary layer on a flat wall in the 
presence of constant mass body force. Let's assumed that the force is directed against 
the flow. The equations describing such a flow are as follows: 


du 3v 
dx dy 


P 




dp dx 
dx dy 


with traditional boundary conditions at the wall and at the core flow. Let us introduce 
an external force potential 


<2>=fx , 


where ,v is a distance between leading edge of the flat and an observation point. 
Then the momentum equation can be rewritten as 


P" 


du 

dx 



d 

- — (/? + <&) 

dx 


dx 

dy 


and it is clearly seen now that the full similarity exists between the positive pressure 
gradient and longitudinal external body force effects. Then the separation criteria (21) 
can be rewritten for new conditions in the form: 
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0 . 


Re e 



( — -/ ) = 0.058. 


In the case with no external pressure gradient the latter can be simplified to 


9 , 

Re 0 '• — — / = 0.058. 

P«r" 

This new separation criterion has been obtained under assumption of the full 
similarity of the pressure gradient and body force effects and. consequently, for 
supersonic flow supposes implicitly that the characteristics length of body force action 
is equal (or the same order of magnitude) as it is for the pressure gradient case, i.e. 
defined by formula (22). From the other hand characteristic length of the body force 
action is determined rather by external reasons such as. for example, a magnetic system 
size. Consequently, the separation criteria should be modified. In order to do this we 
assume here that the separation work is independent on physical reason of separation. 
Then 


rl - 


: dp ' , Pz-Pl , , Pi , N 

: /„ = — ; Ip = P \( — -0 


i dx ' 


l, 


p I 


and using (22) and (23). the modified criteria can be found in the following form: 


' 0.5 


— — t = 0.684 — j Re a : 25 ( M e : - 1 )~ a25 
P c u e~ VP ej 


(24) 


where / is the characteristics length of the region where the body force acts. 

At the point where p—pu the following formula is valid: 

— v = 0. 1 1 (25) 

p, in- 
here pi, — is an ambient atmospheric pressure at altitude h. 

As a confirmation of the criteria obtained above the result of the experimental 
data on boundarv layer separation published in paper [35] can be presented in the form 


ptf 


that coincides practically with (25). 

It was showm in numerical simulation studies [36-38] that the turbulent boundary 
layer separation at the electrode zone outlet happened when the integral MHD 
interaction parameter becomes equal 
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(- 6 ) 


IB < v > 
C < u> 


= 0.043 


where I=j<\>< z> — total channel current. 

<r> — characteristics transversal size of the channel. 
G=<pu><y><:> — characteristics mass flow rate. 

These expressions can be transformed to the form similar to (25) 


— = 0.043 


that is rather closed to the local criteria value. 


Separation by induced current. 

Using (6) the body force needed for separation is defined as 

/ = 0.684) j Re ft - 0 ' 1 25 ( M; - 1 )“°' 25 (27) 

Taking /=0.5m as a minimal linear scale of an external magnetic system one can 
find that 

/= 0.370 MX/m 3 . (28) 

that is very closed to the estimation presented in paper [39]. 

Electromagnetic force is jxB and can be estimated for the induced current case as 
f~ouB l . where o is electrical conductivitv. Thus the needed magnetic induction is 


For rough estimation we use two different semiempirical formula for electrical 
conductivity presented in Figs. 11-12: the first case corresponds to «optimistic» 
evolution of conductivity of PAMIR working fluid seeded with Cs: conductivity was 
calculated with the formula 


' T * p 
a = ct..| — — 

T ; i p„ 


exp(-W(— -— )) 
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where cr x -43.3Sni/m at 7"*=2370K and />*=(). 1 85MPa. JT=2255~K; the second case 

voi i esponds to the COi seeded with C.v used in MHD generator experiments on the 
facility with shock tube [40]: the conductivity was calculated with the formula 

o = ar - v - p " 1 ' : ex p( -IV IT) (31 ) 

where 0 = 0 . /ox 1 0- 4 is the mass fraction of Cs vapor in the working fluid. 

Foi the fust case the magnetic induction level providing the boundary layer 
separation is found as high as 5=8. 4T that seems to be far from the thecnicallv realized 
values for a transportable device. Thus it can be calculated that the controlled 
boundaiv layei separation in the typical nozzle flow is hardlv be provided with the 
induced current only because of relatively low temperature (and conductivitv) in the 
flow. 


Separation by externally applied current. 


From the expression (28) one can define that the value of current densitv is to be 
not lower thany=2.3xl 0 J A/m : for the case of the external magnetic induction is equal to 
2T. 

The power supply needed is estimated with energv equation in the form 


JT y 
P uc ,>— - = — • 

ax a 


or. in an integral form 


l_ : C(T)dT = i ; — — dx. 

' } P"C p 

Using some linearization procedure described in Ref. [30], a total heat power 
delivered in the current affected region is 

O = 2xro J () / " fa d.\ =2rtropMc „ ( T z - T t ) . 

For the minimum current density estimated earlier in this section /=2.3x I OUA/nv 
the total head power is about <9=26MW. This value is approximately equal to 5°» of the 
typical rocket engine power under conditions considered in [30]. 
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CONCLUSION 


( 1 ) The power facility considered with a built-in multipolar MHD generator makes it 
possible to obtain electric power using the MHD conversion of a part of the aas 
flow energy in the peripheral zone of the rocket engine jets. 

(2) The analysis of volumetric MHD effects and output characteristics of a multipolar 
MHDG with electrode modules of different design has revealed that the theoretical 
dependences are in satisfactory agreement with the experimental data obtained in 
the plasma model of a multipolar MHD channel. 

(M An optional design of a power facility has been developed, in which a multipolar 
MHDG is built in the nozzle. The MHD generator utilizes the products of 
combustion ot cryogenic fuel H:;-rO: ; . seeded with cesium. It is demonstrated that 
such a facility can provide the power output of about 10MW. The loss of specific 
thrust of the propulsion unit is relatively low (less than I V). 

(4) The proposed version of a power facility with a built-in multipolar MHDG may 
show promise when the engine utilizes other fuel vapours (e.g.. ox v sen and is 
kerosene seeded with potassium). 

(o) A simplified model and corresponding estimation procedure for MHD separation 
or turbulent boundary layer has been proposed. 

(6) It was round that for typical conditions of rocket engine nozzle flow an induced 
cuii^-nt sepai ation conti ol is hardly reachable due to low boundarv (aver 
temperature and practical limitations for magnetic induction level. 

( . ) An externally applied electrical power can provide the controlling boundary 
separation in a nozzle flow. 
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3. Dependences of the dimensionless averaged quantities <J • E>. <JxB>- of the 
snort-circuit current /y^and the no-load voltage L ;;j on the relative size of the 
electrode moduies in MHD channels with (a) U-shaped electrodes and (b) 
segmented electrodes. 
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Fig-5. Experimental voltage-current U(I) and power .V / characteristic. 
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Fig. 8. Relative friction coefficient against pressure gradient parameter G\ 
(see text for more explanation). 
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Fig.9. Relative value of the parameter n in velocity profile power low against 
pressure gradient parameter Cn (see text for more explanation). 
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Fig. 10. Separation criteria vs Mach number. 
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Fig. 1 1 . Conductivity profile across a boundary layer. 

1 - MHD Generator PAMIR conditions. 

2 - Combustion products of CO-O: Cs seed. 



Fig. 1 2. Effective electromagnetic force profile across a 
boundary layer. 

1 - MHD Generator PAMIR conditions. 

2 - Combustion products of CO-O : Cs seed. 
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